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FOREWORD 


The Carbon Deposition Model for Oxygen-Hydrocarbon Combustion Program (Contract 
NAS 8-34715) was conducted by Aerojet TechSystems (AT) for the NASA-Marshall Space 
Flight Center. The Task VI activities described in this final report were conducted from April 
1989 to May 1990. Fred Braam was the NASA-Marshall Program Manager. At AT, the 
program manager for Task VI was Dr. Sanders D. Rosenberg. The project engineer on Task VI 
was Darby B. Makel. Major contributions to the model developed on the program were made by 
Dr. Ian M. Kennedy of the University of California, Davis (UCD). In addition, injector mixing 
tests were performed at Dr. Kennedy's research laboratory at UCD. Significant contributions 
were made to the program by the following AT personnel: 

Wendelin M. Burkhardt 
Darlene Hunt 
Shirley J. Reed 
Jeffrey C. Schneider 
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ABSTRACT 


The formation and deposition of carbon (soot) has been studied in the Carbon Deposition 
Model for Oxygen-Hydrocarbon Combustion Program (Contract NAS 8-34715.) An empirical, 
1 -dimensional model for predicting soot formation and deposition in L02/hydrocarbon gas 
generators/preburners has been derived. The experimental data required to anchor the model 
have been identified and a test program to obtain the data has been defined. In support of the 
model development, cold flow mixing experiments using a high injection density injector have 
been performed. The purpose of this investigation was to advance the state-of-the-art in 
LC> 2 /hydrocarbon gas generator design by developing a reliable engineering model of gas 
generator operation. The model has been formulated to account for the influences of fluid 
dynamics, chemical kinetics, and gas generator hardware design on soot formation and 
deposition. 
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l.o INTRODUCTION 


This final report describes Aerojet TechSystems' (AT) activities for Task VI: Data Analysis 
and Formulation of an Empirical Model for the Carbon Deposition Model for Oxygen/Hydrocarbon 
Combustion Program, NAS 8-34715. The objective of Task VI was to develop a framework of an 
empirical Low Mixture Ratio Combustion and Carbon Deposition Model and establish a plan for 
validating the model. The model which has been developed is based on phenomenological 
descriptions of soot formation, growth, transport, and deposition. Data to construct the model has 
be derived from three sources. 

1 . Carbon deposition characteristics of L02/methane, LC>2/propane, and LO 2 / RP- 1 
determined from the Task V gas generator tests. 

2 . Open literature; in particular, work conducted during the past eight years. 

3 . Injector mixing data obtained using nonintrusive laser diagnostics on cold flow tests 
conducted with test hardware. 

1.1 BACKGROUND 

Advanced engine studies have indicated that significant mission performance and life 
cycle cost benefits are associated with L02/hydrocarbon propellant engine systems. Potential long 
range applications include a liquid rocket booster for the space shuttle, a block II propellant system 
for OMS and RCS engines on the space shuttle orbitor, an advanced single stage-to-orbit shuttle 
replacement vehicle, and a heavy lift launch vehicle. 

Before the start of the current program in 1982, the development of combustion devices 
for the next generation of booster engines was hampered because of both the lack of critical data 
and the existence of contradictory data for both main chamber combustion devices and prebumer or 
gas generator devices. Low pressure hydrocarbon engine data showed that carbon deposition on 
the chamber wall forms an insulating layer which could be beneficial from a cooling standpoint. 
The extent of carbon deposition on the main chamber walls may determine an engine's operating 
pressure capability by minimizing the combustion heat flux at the chamber wall. 

In prebumers/gas generators, carbon buildup on the turbine blades and nozzles 
increases operation and maintanence costs and may prohibit the use of heavy hydrocarbon fuels for 
reusable applications. The design of LO 2 / hydrocarbon engines requires an understanding of the 
factors controlling soot formation and deposition. However, an incomplete data base on low 
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1.1, Background (cont.) 

mixture ratio combustion characteristics and soot deposition rates for LO 2 /RP-I and L02/methane, 
and the lack of data for L02/propane, have made it difficult to characterize LO 2 /HC combustion for 
the rational selection of the most promising propellant combination and combustor designs for 
future engine development programs. 

In the past, expendable booster engines, operating with LO 2 /RP-I propellants have 
been designed (through costly iterations) to operate with excessive carbon deposition. Future 
reusable engines will not be as tolerant to the performance losses, degraded life, and clean up costs 
associated with soot formation. The first phase of the program, reported in Reference 1 , addressed 
the generation and deposition of carbon using subscale hardware. LO 2 /RP-I was studied at main 
chamber mixture ratios. LO 2 /RP-I, L02/methane, and L02/propane were studied at low mixture 
ratio, gas generator/prebumer conditions. One universal test set-up using the same fine pattern 
triplet injector was used throughout the testing. Carbon deposition during the main chamber 
operation with LO 2 /RP-I was studied for mixture ratios of 2.0 to 4.0 and chamber pressures of 
1000 to 1500 psia (6.89 to 10.34 MPa). Thermal data and visual post-test inspection showed no 
evidence of carbon deposition on the chamber walls. The deposition of carbon on the turbine 
simulator tubes during the prebumer/gas generator testing was evaluated at mixture ratios of 0.20 
to 0.60 and at chamber pressures from 720 to 1650 psia (4.96 to 1 1.38 MPa). The mixture ratios 
covered the bulk combustion temperature range of interest for state-of-the-art turbopump 
machinery, 1300 to 1600 F (978 to 1 144 K). The results show that the carbon deposition rate is a 
strong function of the mixture ratio and a weak function of chamber pressure. Gas generator 
testing with LC>2/propane revealed a threshold mixture ratio for which carbon deposition begins 
and becomes very heavy. Carbon deposition was not detected for LC>2/methane at any mixture 
ratio tested. From the carbon deposition analyses, the turbine drive operating limits were defined 
for each fuel tested. The data indicate that methane is the only hydrocarbon fuel tested that can be 
run without carbon deposition over the desired gas generator operating temperature range. 

The results from the first phase of testing for the carbon deposition and gas tempera- 
tures as a function of mixture ratio are summarized in Figure 1. The curves on each plot indicate 
the measured gas temperature as a function of mixture ratio tested for each fuel. Superimposed on 
each plot is the desired temperature range for the operation of state-of-the-art turbine drives. The 
highlighted area indicates the region where operation for each fuel will not incur or at least mini- 
mally incur carbon buildup. The intersection of the highlighted area with the area delineating the 
desired temperature range indicates the region of acceptable performance for a gas generator for 
each fuel. Figure 1 indicates that LO 2 /RP-I cannot be operated in the desirable temperature 
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Figure 1. Gas Generator Operating Limits Defined in First Phase of Program 


1.1, Background (cont.) 

range for gas generators without incurring undesirable carbon buildup. LC> 2 /propane can be 
operated in 'the desired temperature range up to a maximum of 1500 F (1088 K). Operation with 
L02/methane is unrestricted over the desired gas generator operating temperature range. 

The second phase of the program, reported in Reference 2, extended the L02/methane 
database by investigating fuel purity and injection density effects on carbon deposition at gas gen- 
erator conditions. To test the effect of fuel purity, L 02 /Hquidified natural gas (LNG) was tested 
using the same fine pattern injector as the earlier testing. The LO 2 /LNG gas generator/prebumer 
testing was performed at mixture ratios between 0.24 and 0.58 and chamber pressures from 840 to 
1 370 psia (5.8 to 9.4 MPa). A course element triplet injector was constructed to permit testing at 
approximately full scale injection densities (10 times the injection densities used in the earlier 
testing). The new injector was used for high injection density LC>2/methane testing at mixture 
ratios from 0.23 to 0.81 and chamber pressures from 925 to 2210 psia (6.4 to 15.2 MPa). No 
detectable carbon deposition was measured for for the LO 2 /LNG and the high injection density 
L(> 2 /me thane tests. 

While the results of the earlier phases of the program successfully characterized the 
carbon deposition differences among the three hydrocarbon fuels tested, sufficient data to develop 
a model to predict soot formation and deposition were not obtained. It was recommended at the 
end of the second phase of testing that the effects of hardware geometry, injector design, and 
combustion kinetics be investigated as potential means of reducing the high soot deposition rates 
measured for LO 2 /RP-I and L02/propane. The objective of the investigation would be to develop 
an engineering model for gas generator design which can define conditions leading to minimum 
carbon deposition. 

Several advances in the application of laser diagnostic techniques to practical combus- 
tion devices have occurred in the past 10 years. Application of these techniques to gas generator 
combustion should enable the combustion and carbon deposition to be empirically modeled using 
detailed measurements of local flow and combustion parameters. Pnor attempts to model gas gen- 
erator soot formation and deposition have proven inadequate for design purposes because the soot 
concentration in the flow, the most important parameter for the model could not be measured. The 
approach taken in this study has been to formulate a model relating soot formation and deposition 
to flow and combustion parameters such as propellant type, local wall shear stress, turbulence, 
velocity gradients at the interface of propellant streams from individual fuel and oxidizer injector 
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1.1, Background (cont.) 

elements, and temperature gradients. This approach should yield a model that predicts soot 
formation and deposition in terms of parameters which can be controlled through hardware design. 

1.2 OBJECTIVES 

The first objective of the program was to formulate a Low Mixture Ratio Combustion 
and Carbon Deposition Model based on available data and realistic models for soot formation, 
growth, transport, and deposition. The second objective was to define an experimental program, 
based on the use of nonintrusive laser diagnostics to provide the key information needed for verifi- 
cation of the model. The model is intended to be a tool for the design of gas generator components 
such as injectors and combustion chambers, based on a minimum soot formation criteria. The 
model will also predict the total soot yield and local carbon deposition as a function of operating 
conditions and geometry. The experimental program is based on a progression from laboratory 
scale experiments to high-pressure, single-element injector experiments, and eventually to subscale 
experiments using the existing test hardware from NAS 8-34715. 
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2.0 SUMMARY 


An empirical model of carbon formation and deposition applicable to LC>2/hydrocarbon gas 
generators/prebumers has been formulated. A literature review has been conducted to review the 
state-of-the-art in soot formation and deposition modeling. Previous models and test data were 
evaluated and have been used to construct the model. The purpose of this modeling effort has been 
to develop a unified methodology for predicting soot formation and deposition rates when 
designing L02/hydrocarbon gas generators. The equations derived for the model can be pro- 
grammed into a computer code. The empirical constants in the model have been defined and a test 
program for quantifying these constants and validating the model has been prepared. 

The most important physical mechanisms controlling soot formation and deposition at low 
mixture ratio gas generator conditions are contained in the model. These mechanisms are: ( 1 ) 
inception of solid soot particles from gaseous combustion products, (2) increase in the size and 
mass of soot particles due to surface reactions with gas phase hydrocarbons and collision with 
other soot particles, (3) deposition via thermophoresis (i.e. "mass transfer" down a temperature 
gradient,) and (4) reentrainment and shearing of previously deposited soot at high flow velocities. 
The derived equations for the model are based on a one-dimensional representation of the flow in 
the gas generator and turbine drive assembly. The concentration of soot (expressed as the volume 
of soot per unit volume of flow) is calculated as a function of time. A chemical kinetics model (yet 
to be defined) for each propellant combination is used to calculate the axial temperature profile and 
concentration of unbumed C2 and C3 hydrocarbons at stations down the length of the chamber. 
Once the axial profile of soot concentration and temperature are known, the deposition rate of soot 
on the wall at each station is calculated. The deposition rate is determined using a mass transfer 
analogy to the local heat transfer. The local wall heat transfer coefficients are calculated using 
existing correlations and are then transformed into soot mass transfer coefficients, which account 
for the effects of thermophoresis and reentrainment. 

A major unknown issue in the modeling of soot formation at low mixture ratio combustion 
conditions is the effect of injector mixing. The mixing characteristics of the injector determine 
whether the combustion is primarily premixed (kinetically limited) or diffusion dominated (mixing 
limited). Depending on the injection pattern, a gas generator may exhibit varying degrees of each 
type of combustion. The chenical reactions leading to the formation and growth of soot particles 
are slow reactions relative to the gas phase combustion reactions. These reaction occur on time 
scales approximately equal to those of the turbulent mixing, and therefore are profoundly affected 
by the fluid dynamics in the combustion zone. Cold flow tests to assess the mixing characteristics 
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2.0, Summary (cont.) 


of the high injection density injector (considered a conventional gas generator injector design) have 
been conducted. The tests used molecular scattering of laser light to visualize 2-D cross sections 
of the flow near the injector face The tests reveal large zones of unmixed fuel at more than 20 
injection orifice diameters downstream of the injector face. Based on these results a particle 
nucleation rate term dependent on the mean axial turbulent mixing rate has been included in the 
model. 
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3.0 CONCLUSIONS AND R ECOMMENDATIONS 


3.1 CONCLUSIONS 

Currently, there is no reliable methodology for predicting soot formation and deposition 
in LC> 2 /hydrocarbon gas generators. The processes governing soot formation and deposition are 
highly dependent on hardware configuration (injector pattern, chamber contour) and operating 
conditions (MR and Pc.) Therefore, existing test data cannot be directly extended to new designs. 
The viability of heavy hydrocarbon fuels (RP-1 and propane) for future missions is limited by the 
present lack of knowledge as to how to design gas generators which produce minimum soot 
deposition. The initial step in filling this technology void has been taken on this program with the 
development of the framework of an engineering model for gas generator design. 

The equations for an empirical Low Mixture Ratio Combustion and Carbon Deposition 
Model have been formulated. The model represents a significant advancement over previous 
approaches to predicting the soot formation and deposition tendencies of L02/hydrocarbon gas 
generators. The model accounts for the fundamental physical processes controlling soot formation 
and deposition. The advantage of this type of model is that it is applicable to different hardware 
designs and should be valid over a wide range of operating conditions. The major tradeoff for this 
type of model is that it requires data types beyond what has been collected to date in order to quan- 
tify parameters and validate the model. Specifically, the soot particle inception rate as a function of 
injector mixing, the soot surface growth rate for the candidate hydrocarbon fuels at elevated pres- 
sures, and the soot layer adhesion strength must be determined. These data can be obtained using 
proven laser-based combustion diagnostic techniques on subscale tests such as those conducted 
during Task V. 

In addition to the development of the model to predict soot formation and deposition, 
several hypotheses regarding how to reduce or entirely eliminate soot deposition have been 
deduced during the program. These hypotheses are: (1) the suppression of soot particle nucleation 
through the use of very rapid mixing of fuel and oxidizer, (2) the suppression of soot surface 
growth reactions through the use of fuel additives, (3) the reduction of deposition through tailored 
design of the wall temperature profile. Exploiting one or a combination of these mechanisms 
through innovative hardware design may result in a "soot free" gas generator design. 
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3.0, Conclusions and Recommendations (cont.) 

3.2 RECOMMENDATIONS 

The program has provided considerable data and the framework of an engineering 
model for the prediction of soot formation and deposition. A systematic program should be initi- 
ated to continue the development of the model and perform the required model validation tests as 
described in the technical summary (Section 4.0). The empirical inputs required by the model 
should be determined through a progression of controlled experiments. The major component of 
the test program should be direct measurements of the soot concentration and deposition using the 
existing Task V test hardware and nonintrusive laser diagnostics. The modular design and known 
combustion stability characteristics of the existing program hardware make it ideal for validation 
testing. Laser diagnostic access for the required measurements can be readily accommodated by 
means of special "diagnostic flanges" which could easily be incorporated into the hardware test 
configuration. Following the model validation tests, the model should be used to design a high 
pressure L02/hydrocarbon gas generator. This gas generator would be designed to minimize soot 
formation and deposition and to serve as a baseline design for a future LCfyhydrocarbon engine 
development program. 

In parallel with the model validation testing, the existing 1-D model should be extended 
to 2-D or 3-D as a submodel to be integrated into a computation fluid dynamics (CFD) code. CFD 
is emerging as a necessary design tool for high pressure, high performance engine design and 
analysis (e.g. SSME product improvement and ALS.) A CFD-based soot formation and 
deposition model is potentially the ultimate analysis tool for detailed gas generator design. Such a 
model should greatly reduce the inherent limitations of an empirically based 1-D model. CFD 
based models of soot formation and transport in turbulent flames have proven successful and 
appear readily adaptable to current rocket engine combustion codes. Concurrent development of a 
CFD based model and the empirical model validation tests would maximize test effectiveness and 
reduce the overall development time of a useful CFD model. 
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4.0 TECHNICAL DISCUSSION 


4.1 FUNDAMENTALS OF SOOT FORMATION AND DEPOSITION 

Solid carbon (soot) deposits found in gas generators can be formed by two mecha- 
nisms. The first mechanism is pyrolysis of liquid fuel due to droplet heating during combustion. 
This process, known as coking, results in solid particulates ranging in diameter from 50 to 200 
microns. This is the same phenomena which causes carbon deposits when hydrocarbon fuels are 
used to regeneratively cooled or film cooled thrust chambers. Direct coking of the liquid fuel 
provides a small contribution to the total solid carbon formation under gas generator combustion 
conditions. The second, and most important mechanism, is soot formation through gas phase 
reactions. During this process, nucleation of small soot particles occurs in the fuel rich side of the 
combustion zones. These particles are very small in diameter, typically 1 to 10 nm, but rapidly 
grow in size to about 200 nm as they react with excess hydrocarbon molecules. The surface 
growth reactions typically account for 90 percent of the solid carbon mass deposited in combustors 

[3]. 

The soot formation process in a gas generator is illustrated in Figure 2. Discrete injec- 
tion locations for fuel and oxidizer create hot, fuel rich zones comprised of gaseous fuel and partial 
vaporized droplets near the local flame fronts where the mixture ratio is stoichiometric. At low 
mixture ratio conditions, the fuel rich zones must exist. An injector with high mixing rates will be 
characterized by flame fronts with high surface area and the fuel rich volume of the flow distributed 
into many small zones. An injector with low mixing is characterized by a flame front with rela- 
tively low surface area and large zones of unburned fuel. The unbumed, vaporized fuel molecules 
undergo thermal decomposition resulting in the formation of polycyclic aromatic hydrocarbons 
(PAHs.) The PAHs are known to be precursors to soot formation although the detailed kinetics of 
the formation process are unclear [3]. Soot particles are nucleated within the hottest portions of the 
fuel rich zones. The residence time of a particle in a hot zone is significantly less when the mixing 
is very rapid. The particle will continue to grow while it is in a region of the flow where the local 
temperature exceeds 2500 F (1558 K.) Once the soot particle is beyond the flame zone and its 
temperature is reduced to the typical gas generator bulk exit temperature of 1300 to 1600 F (926 to 
1084 K), the particle stops growing. This effect is explained by the fact that the surface growth 
reactions between the gas phase hydrocarbons and the soot particles, which result in the addition of 
solid carbon to the particle, are highly temperature dependent. The temperature dependance of the 
reactions explains why at extremely low mixtures (less than 0.25) little soot is formed. At these 
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Figure 2. Soot Formation Processes 


4.1, Fundamentals of Soot Formation and Deposition (cont.) 

conditions, the combustion is so oxidizer deficient that the flame regions are necessarily very small 
and the heat release is too low to accommodate soot particle nucleation or growth. 

The soot deposition processes within gas generators and turbomachinery components 
are illustrated in Figure 3. Once the combustion gases are thoroughly mixed (usually by a turbu- 
lence ring) the flow can be characterized as a hot, nonreacting, particle laden flow. The soot 
particles in the flow are approximately 200 nm in diameter [3]. Particles of this size are too large to 
be significantly influenced by diffusion (concentration gradients) and are too small to have suffi- 
cient inertia to deviate from the fluid streamlines. The process primarily responsible for soot 
deposition has been shown to be thermophoresis ( mass transfer down a temperature gradient ) 

14]. Thermophoresis arises from the fact that a particle in a flow with a temperature gradient is 
struck with high energy (i.e. high velocity gas molecules) on the hotter side of the particle than on 
the colder side. The net effect of the imbalance of forces on the particle is a net force propelling the 
particle across streamlines from a hot gas toward a colder surface [5]. Even when the combustion 
gases and wall temperatures vary by as little as 5 %, thermophoretic deposition is the dominate 
deposition mechanism, over a factor of 10 greater than diffusional deposition. 

As illustrated in Figure 3, the soot deposition patterns in a gas generator show the 
heaviest deposition in areas of high heat transfer such as stagnation surfaces and nozzles (high 
velocity regions). At sufficiently high velocities, however, such as in turbine nozzles, the fluid 
wall shear stress can exceed the adhesion strength of the soot particle to the wall. When the soot 
adhesion strength is exceeded, reentrainment of soot particles will occur and the net deposition will 
be reduced in the region. Cyclic differential pressure fluctuations, such as those measured during 
testing of the Titan I gas generator [1], were probably due to an increasing wall shear stress as the 
deposition layer built up and the local flow velocity increased, due to a reduction in the flow area. 
The soot layer would then build up again once some of the layer was sheared away and the local 
velocity was reduced so that the wall shear stress was less than the average soot layer adhesion 
strength. The important consequence of the dependance of deposition on wall shear stress is that 
different gas generator flow geometries will experience different "net deposition," even when the 
amount of soot formed is the same. 
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4.2 MODEL FORMULATION 


Approaches for modeling soot formation and deposition are separated into three cat- 
egories representing different levels of sophistication and applicability. The first level (Level I) 
consists of single parameter correlations of experimental data. The second level (Level II) consists 
of one dimensional, lumped-parameter methods in which the flow equations, combustion kinetics, 
soot formation and deposition are determined at discrete locations along the length of the gas 
generator and turbine inlet. The most sophisticated modeling (Level III) involves solution of 2-D 
or 3-D conservation equations for momentum, energy, and soot volume fraction (i.e. volume 
occupied by soot particle per unit volume of space in the flow field.) 

The Level I type models represent existing design tools and describe the complicated 
physics occurring within a gas generator inadequately [1]. The Level II type models are 
empirically based analytical models which are useful for design provided sufficient data is obtained 
to properly adjust the empirical constants. The Level III type models represent the ultimate design 
and analysis tool, however such modeling at present is limited by the uncertainty of CFD solutions 
of the two-phase turbulent combustion typical of gas generators. The Level III type models should 
be useful in modeling controlled tests which may be required to completely develop a Level II type 
model. A summary of the various models and their relative advantages and disadvantages are 
summarized in Figure 4 and discussed below. 

Level I - Level I type models were extensively summarized in the first Interim 
Summary Report [1]. These models include the Seader and Wagner carbon deposit correlation and 
the Rocketdyne correlation of carbon deposit thermal resistance. These types of correlations are 
based on single parameters (such as mixture ratio and mass flux of combustion gases) and have 
been developed based on limited data sets. These models are useful from the standpoint of 
indicating trends in soot formation (e.g. gross MR dependence), but provide no reliable data for 
developing new gas generator designs. Level I type models offer little potential for advancing the 
state-of-the-art in gas generator design. 

Level II - Level II type models appear to have the greatest potential as a near-term 
engineering models for soot formation and deposition. Farmer et al [6] and Jensen [7] developed 
Level II type models for predicting soot formation and total soot yield from well mixed jet stirred 
reactors and methane combustion, respectively. The soot nucleation, growth, and oxidation rates 
were expressed in the form of Arrhenius type rate equations. The effect of local mixing was not 
included in these models. However, the soot residence time in the combustion chamber was 
included. The empirical constants were adjusted so that an acceptable fit to available experimental 
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Figure 4. Modeling Approaches for Low Mixture Ratio Combustion and Carbon Deposition Model 















4.2, Model Formulation (cont.) 


data was achieved. The basic form of the soot formation and growth equations provides a basis 
for the development of a Level II Low Mixture Ratio Combustion and Carbon Deposition Model. 

The prediction of soot deposition is analogous to chamber heat transfer for which 
lumped parameter, 1-D models have been proven as adequate design tools once a sufficient amount 
of experimental verification of empirical constants has occurred. An example of a Level II heat 
transfer model is the NASA LeRC HOCOOL program for calculating thrust chamber heat transfer 
with regenerative and film cooling [8]. Correlations for mass transfer coefficients based on local 
flow properties including the effect of soot particle thermophoresis have been developed by 
Gokoglu and Rosner [9] and Rosner [10] to predict particle deposition in gas turbines. The net 
deposition of soot to the wall is ultimately controlled by the adhesion strength of the soot layer 
versus the shearing action to the high speed flow near the wall. The empirical adhesion strength 
model for small particles, developed by Zimon [11] and Barengoltz [12] (see Section 4.2.4 for 
details), is applicable to soot deposits provided the mean adhesion force is known. 

Aerojet's Fuel Rich Combustion Model (FRCM) [13] is an example of a Level II type 
model using a chemical kinetics scheme to calculate the average temperature for fuel rich LOX/RP- 
1 combustion. This model does not predict soot formation, deposition, or include the effects of 
mixing processes and chamber geometry. The reaction scheme used in the FRCM could easily be 
adapted for use in a Level II Low Mixture Ratio Combustion and Carbon Deposition model in 
order to account for the pyrolysis reactions of LOX/RP-1. 

Level III - The Level III models represent the highest level of sophistication and ulti- 
mately the most powerful modeling capability. The objective of such models is to predict the 
formation, growth, and transport of soot throughout the flow field. Detailed kinetic models of soot 
formation have been studied by several researchers such as Frenklach et al [14]. In general, such 
reaction schemes (involving on the order of 500 to 1000 equations) have been unsuccessful in 
predicting soot formation and are computational intensive. The current trend in Level III type 
modeling uses empirically derived relationships between the local mixture ratio, flow field, and 
temperature field to the amount of soot formed. An early attempt at this type of modeling was 
performed by Gore and Faeth [15] for ethylene diffusion flames using an empirical correlation 
between the soot volume fraction and the local mixture ratio. This correlation was obtained from 
measurements made in a laminar axisymmetric diffusion flame. The problem with this correlation 
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4.2, Model Formulation (cont.) 

is that it ignores the effect of the soot particle history and residence time in the flame. It also 
ignores the effect of mixing rates on the soot chemistry. 

Magnussen and Hjertager [16] and Magnussen et al [17] incorporated a model for soot 
formation into an eddy dissipation model for turbulent combustion. They used the kinetic scheme 
of Tesner et al [18] for soot formation from air/acetylene combustion by calculating the mass of the 
soot and mean particle number density. The empirical constants were adjusted to give good 
agreement between the computations and their experiments. 

Recently Moss et al [19] proposed a model which incorporates the essential physics 
i.e., a nucleation rate, a surface growth rate and a bum-out or oxidation rate. Their rate constants 
were not known and were determined by fitting their model to measurements in a laminar diffusion 
flame. The model which they propose may be somewhat difficult to implement in to turbulent flow 
calculations because of the usual closure problems for the chemical source terms in conventional 
models of reacting flows. Methods tusing transport equations for probability density functions 
should overcome the turbulent closure problem. 

Kennedy et al [20] have recently developed an improved Level III type model which 
can be easily incorporated into existing turbulent combustion codes. It should be useful for 
modeling specific tests to define the empirical constants for a Level II type model. In the absence 
of adequately reduced chemical schemes for soot formation, the approach is to use the available 
empirical information which laser diagnostic measurements have provided. In addition, this model 
is based on a satisfactorily accurate numerical scheme which does not add substantially to the 
already significant computational burden that is faced in predicting practical turbulent flow fields. 

In Kennedy’s model, the local mixture ratio is the primary quantity which is calculated. 
The temperature, density, and the gas composition are determined as functions of the local mixture 
ratio. The local mixture ratio is also used to determine the soot volume fraction indirectly. The 
soot volume fraction is not a function of the local mixture ratio, but rather the rates of nucleation, 
surface growth, and oxidation are functions of local mixture ratio. The conservation equation for 
the soot volume fraction, <|> , is then 


dp / \dp 1 
pi, — +p(v + v T )— = - 


dr 


f dp 

rpD - 9 


V 


s dr 


\ + pw K +pw g - pw 0 


(1) 
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4.2, Model Formulation (cont.) 

where w ,w , and w q are the rates of soot volume formed by nucleation, surface growth, and 
removedby oxidation respectively. A thermophoretic radial velocity, v T of soot particles is 
included and is very important in determining soot deposition rates. This model has been used 
successfully to predict soot formation and transport in laminar and turbulent flames. 


4,2.1 Qvprall Methodology 

The model chosen for development in this program was the near-term, tractable 
model (i.e. Level II model) consisting of an empirically based 1-dimensional representation of the 
gas generator combustion and flow. The essential elements of the model are (1) combustion kinet- 
ics submodel, (2) soot formation submodel, (3) 1-D flow model, and (4) a soot formation model. 
The outputs of the model are axial profiles of temperature, gas composition, soot yield, pressure, 
wall shear stress, wall thermal resistance, deposit thickness, and deposition rate. The elements of 
the model are illustrated Figure 5. The basic inputs to the model are propellant type, a mixing 
parameter characterizing the injector, flowrate, mixture ratio, chamber pressure, and chamber 
geometry. From these inputs, a soot formation model calculates the soot particle nucleation and 
growth rates. The concentration of combustion products and bulk temperature are provided by a 
kinetic model consisting of a reduced set of gas phase combustion and pyrolysis reactions The 
soot formation model is iterated with a 1 -D flow model to determine the amount of soot (soot vol- 
ume fraction), bulk temperature, and wall shear stress down the length of the chamber. 

In this model, the soot formation process is decoupled from the deposition 
process, except to the extent that soot build up over time may cause an increase in chamber pres- 
sure. Therefore, once the flow conditions and soot volume fraction profile have been specified 
down the length of the chamber, the deposition model is used to predict the thickness of the soot 
deposit. The input to this part of the model is the wall shear stress, soot volume fraction, and wall 
temperature at specified stations down the length of the gas generator. The soot deposition at each 
station is calculated in two steps. First, a mass transfer coefficient for soot is used to calculated 
the maximum possible deposition. The net deposition is then determined by calculating a "sticking 
fraction" of soot on the wall by comparing the soot layer adhesion strength to the wall shear 
stress. 
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Figure 5. Elements of the Low Mixture Ratio Combustion and Carbon Deposition Model 















4.2, Model Formulation (cont.) 


Models currently exist for the chemical kinetics and 1-D flow which are 
directly applicable with only slight modification. The kinetics submodel is used to predict only 
combustion temperature, but not soot nucleation. Therefore, the kinetics model can be replaced 
with empirical data at this level of modeling. The emphasis during this program was placed on 
formulating the basic equations of the soot formation and deposition models. 

4.2.2 Identificatio n of Model Parameters 

From the results of the literature review key parameters affecting soot forma- 
tion and deposition at low mixture ratio conditions have been identified. An initial screening of the 
effects of combustion parameters on soot formation and deposition was conducted in the first 
phase of the program and was presented in Reference [1]. That database has been extended to 
include candidate parameters for the model . The results of this review are the parameters which 
have been incorporated into the models for combustion, soot formation, and deposition. A 
summary of key parameters is shown in Figure 6 where the effect on soot formation and 
deposition is indicated. 

4.2.3 Soot Formation Model 

A soot formation model has been derived which accounts for the contribution 
of the inception rate of small (1-10 nm diameter) soot particle nuclei and growth/coagulation of the 
particles as they traverse the gas generator. Because the soot particle inception rate is one of the 
slowest reactions relative to gas phase combustion reactions, it is very sensitive to fluid mechanics 
(i.e. mixing characteristics) of the fuel and oxidizer. Figure 7 illustrates the potential effect of 
mixing (in terms of a relative velocity gradient between spray fans) on particle formation rate. The 
data shown on the figure were obtained in a counter flow diffusion flame where the mixing rate 
was precisely controlled [24]. At sufficiently high mixing rates, there is insufficient time for soot 
particles to form (nucleate) and the overall soot formation is suppressed. The soot particle number 
density at inception drops from a typical value of lO 1 ^ particles/cm 3 very rapidly to nearly zero. 

At low mixing rates, well below that sufficient to cause soot suppression, the dominate processes 
affecting soot formation are the soot surface growth reactions and residence time of soot particles 
in hot regions of the flow. Under these conditions, total soot yield is relatively insensitive to the 
precise number density (i.e. variations of a factor to 10 to 100) at formation because the rapid 
coagulation and growth of "young" soot particles controls the surface area available for the addi- 
tion of solid mass to the particle. The empirical model which has been derived to account for these 
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Figure 6. Parameters Affecting Soot Formation and Deposition 
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4.2, Model Formulation (cont.) 

processes uses measured values for the particle inception rate and the soot surface growth rate for 
particular hydrocarbon fuels. The derivation of the model is described below . 

First, the effect of soot surface growth reaction (i.e. the conversion of gas 
phase hydrocarbons to soot) is derived. The soot volume fraction is defined as the volume of soot 
per unit volume of flow as given by eq.2. 


(2) 

N= Soot particle Number Density {particles/cm 3 ) 
d = Soot Particle Diameter 

P 

The effect of the surface growth reactions is to increase the mass and therefore the volume of soot. 
The change in mass as a function of time is therefore: 

—j~ = Njid 2 P Cp s (3) 

dt 


C= Soot Surface Growth Rate 
p = density of soot particle=1.9 gm/cm 3 


The soot surface growth rate, C, represents an average of all the reactions causing growth of soot 
as shown by eqs. 4 and 5. 


(C n H m ) 4 -^C, 


if— Y exp j 

E a 

1 P 


\o) 

R 

, T « T 'j 



(4) 

(5) 


where: a = rate constant 

c/o= local carbon to oxygen ratio 
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4.2, Model Formulation (cont.) 


E = activation energy 
T c = bulk gas temperature 
T re{ = reference temperature 


An overall rate is used in the model because of a lack of detailed kinetics models. This overall 
reaction rate can be directly measured with the fuels of interest in controlled combustion experi- 
ments. The soot surface growth rates for hydrocarbon/air diffusion and premixed flames have 
been measured and for a given fuel, a single step rate equation such as eq. 5 has proven valid over 
a wide range of mixture ratios [27]. 

Substituting for dp in eq. 3 with eq. 2 and dividing by p,, the rate of change of 
the soot volume fraction as a function of time is obtained in eq. 6. 

£-(«!»' W*C (6) 


Next, the effect of soot particle coagulation which changes the soot particle number density in time, 
N(t), must be derived. Following the derivation of Hinds [28], first the diffusion of particles to a 
single particle is considered and then extended to all the particles in the unit volume. For an 
individual particle the rate of collisions per unit time per unit area is governed by Fick’s law of 
diffusion (eq. 7.) 


J = 


dN 

D— 

dx 


(7) 


D= Diffusion Coefficient 


J-Flux of particles to the surface 


The collision surface is assumed to be twice the soot particle diameter. The rate of collisions 
between the selected particles and the surrounding particles as a function of time is equal to the 
product of the collision surface area and the flux of particles. 


dn a .2 ^ dN 
— = 4 ndp D — 
dt dx 


( 8 ) 


n = number of collisions 
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4.2, Model Formulation (cont.) 


The concentration gradient at the surface is obtained by solving eq. 7. 


dN _ -2 N 

dx d p (9) 


The rate of collisions for all the particles is equal to N/2 times the rate of a single particle. (The 
factor of 2 prevents double counting of collisions.) The rate of collisions per unit volume is: 


^ = 4 n d DN 2 
dt 


00 ) 


n c = collisions per unit volume 

Assuming that for each collision of a soot particle there is a reduction of one in 
the number of particles in a unit volume, the number concentration is reduced by one. The change 
in number concentration is exactly equal to the collision rate as shown in eq. 11. 

dN _ dn c 
= 


For the size range of soot particles, 10 to 200 nm diameter, the product dp is nearly constant, 
therefore, the time rate of change of the number concentration is eq. 12. 



( 12 ) 


where j5=4nd p D 


The term, f), is called the particle collision coefficient. The particle number density as a function of 
time is obtained by integrating eq. 12, assuming an initial particle number density at formation of 

N a . 


N«) = 


i+PKt 


(13) 


An expression for the soot volume fraction as a function of time (or location in 
the gas generator) can be derived by substituting eq. 13 into eq. 6 and integrating. The resulting 
expression is: 
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4.2, Model Formulation (cont.) 


m 



c 



{(A£ fit+ 1) 2 ' 3 - l}l 


(14) 


Equation 14 provides a means of calculating the axial distribution of soot in the gas generator. The 
primary variables in the model are the number density of particles at inception, N (from the parti- 

o 

cle inception rate), the soot surface growth rate, C, and the particle collision rate, (5. The first two 
of these variables are quantities which will eventually have to be anchored experimentally. The 
collision rate of soot particles can be calculated from existing theory and is essentially independent 
of fuel type. The effect of pressure and temperature are contained in the soot surface growth rate 
term. As the flow mixes and the the temperature of the gas becomes uniform the surface growth 
reactions stop and the soot volume fraction becomes a constant. The effect of mixing is contained 
in the number density at particle inception, N , as discussed above. 

o 

For the one dimensional streamtube model approach, all parameters are repre- 
sented by bulk, average values. If the propellants are injected premixed, then the one-dimensional 
model provides a very good approximation for the flow. Discrete and separate injection of the fuel 
and oxidizer adds two features, mixing and local temperature variations, which are very important 
for predicting soot formation. The first feature, mixing, affects the soot particle inception rate. The 
second feature, local temperature varations in the flow, provide high tempeature regions where 
soot particles can grow. Each of these features impacts the total amount of fuel converted solid 
carbon. In order to account for these processes in the model, empirical correlations for the panicle 
number density at nucleation, N 0i and the "effective" combustion temperature used in eq. 5, based 
on the injector mixing efficiency, £ m , must be developed. The basis for the correlation is the 
assumption that a high mixing efficiency that is achieved in a short time (or short chamber length) 
implies rapid fuel and oxidizer mixing and the "effective" combustion temperature must quickly 
approach the bulk temperature. Sufficient data does not exist to quantitatively define this 
correlation, however, the proposed form of this correlation is given in eq. 15. 

Ar„ = £ . W a.Ar 1 .ex P (-c,£) (15) 

where: £jx) = axial mixing efficiency profile 

v= average velocity 
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4.2, Model Formulation (cont.) 


L= chamber length 

The "effective" combustion temperature is then related to the bulk temperature by a correlation ol 
the form in equation 16. 


T c — T hM * 


A 




(16) 


where: T = Average Temperature from kinetic model 


The mixing efficiency can be input from experimental data or calculated for 
conventional injector element configurations using the standard JANNAF performance prediction 
methodology [29]. The empirical constants appearing in eqs. 17 and 18 ,Kj t K 2 , and cj , can be 
anchored for different injector designs by measurements for the soot particle size and soot volume 
fraction at several axial positions in a prototype gas generator. Details of the measurement 
techniques and test requirements are described in Section 4.4. 

The model has been used with the 1986 LC>2/propane tests as test cases. Table 
I contains a summary of the L02/propane Task V data together with the calculation of the soot vol- 
ume fraction (<|>) for each test. The calculation was made assuming a soot surface growth rate of 

C = 5.0e -07(C / 0) 2 e * T « T ‘ 


= 20 


kcal 

mole 


T rtJ = 1700K 

which is typical for premixed hydrocarbon/air combustion systems with C 2 or C 3 type fuels. A 
average particle collision rate was calculated for 0.1 micron soot particles at MR= 0.35, 1000 psia 
and Tc=1400 F. The purpose of these calculations is to see if the proposed model exhibits the 
same trends as those observed in the testing. There is no way to infer the value of the soot volume 
fraction directly from the Task V data. 

Figure 8 shows the calculated soot volume fraction as a function of mixture 
ratio. In general the soot volume fraction in the range of interest (M.R. 0.3 to 0.5) is very sensi- 
tive to mixture ratio. Scatter in the data is due to the fact that the final soot volume fraction is 
dependant on residence time in the gas generator. The absolute value of the soot volume fraction 
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TABLE I. Calculation Of Soot Volume Fraction For Task V L0 2 /Propane Tests 
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4.2, Model Formulation (cont.) 

appears to be of the correct order of magnitude. Sooting flames typically exhibit soot volume frac- 
tions ranging from 1.0e-07 to 1 .0e-05. Inspection of Figure 8 shows over an order of magnitude 
increase in the soot volume fraction between MR=0.30 and MR=0.40. This increase correlates to 
the measured increase in pressure drop across the turbine simulator during testing as shown in 
Figure 9. The soot deposition rate is directly proportional to the soot volume fraction, therefore a 
large increase in soot production is expected to produce large change in the indicated deposition. 

4.2.4 Soot Deposition Model 

The net soot deposition is the difference between the gross deposition rate and 
the removal (reentrainment) rate. To obtain the gross deposition rate, anchored heat transfer 
coefficients are used to calculate mass transfer coefficients. The mass transfer coefficients are then 
modified to account for thermophoresis ("mass transfer down a temperature gradient.") The effect 
of thermophoresis on soot deposition is similar to that of wall suction on heat transfer. For soot 
particles which range in size from 10 to 200 nm, deposition via thermophoresis is several orders of 
magnitude greater than diffusion or inertial impaction. The removal rate is determined separately 
from the flow conditions at the wall. The soot removal rate is determined from a balance of the 
strength of the soot particle adhesion to the wall versus the fluid shear stress at the wall. 

The gross thermophoretic deposition is calculated using the correlation pro- 
posed by Rosner et al [4,20]. This correlation is based on the assumption that the thermophoreti- 
cally enhanced deposition of soot is proportional to the “normal” (i.e. diffusion only) mass transfer 
coefficients, as shown in eq. 17, 

(17) 

where the normal Stanton number for mass transfer is defined in eq. 18. 


The first function in eq. 17 accounts for the transport of soot to the wall due to the overall 
temperature gradient, which is analogous to “wall suction” on heat transfer. The second function, 
f 2 (D), accounts for the fact that the temperature gradient is not linear, and has a nonzero second 
derivative. The effect of a nonlinear second derivative is that the wall acts like a sink for soot par- 
ticles due to a thermophoretic acceleration of particles towards the wall. Based on laminar deposi- 
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4.2, Model Formulation (cont.) 


tion measurements and correlation with turbulent flow CFD calculations the functional form of the 
relations has been found by Rosner [4] to be as shown in eqs. 19 and 20. 


/,(*) = 



(19) 


f 2 (D) = e 


-D 


The variables B and D are defined by eqs. 21 and 22, 


B = -(0.55Le) w |^ 


T.-T w 

\ T w J 


D = (0.55Le) w 


( T _ T ^ 

1 e 1 w 


( 20 ) 


( 21 ) 

( 22 ) 


where: T, = bulk gas temperature 

T w = local wall temperature 


where, the Stanton number for heat transfer, St h , and the Lewis number, Le, are defined based on 
bulk property data as shown in eqs. 23 and 24. 

Nu 


St„=- 


RePr. 


hd 

Nu- — 

k 


(23) 


where: 


Re = 



Pr = 


a£t. 

k 


h= local heat transfer coefficient 
k= bulk gas thermal conductivity 
C p = bulk gas specific heat 
)i= bulk gas viscosity 



( 24 ) 


RPT/F.0068 .66 


32 


4.2, Model Formulation (cont.) 


The local Nusselt number, Nu, can be readily obtained from existing heat transfer correlations or 
derived from combustion chamber heat flux data such as that obtained during the Task V testing. 
Finally the gross deposition rate can be expressed as shown in eq. 25.: 


J deposit QtP V ^m 



(25) 


The deposition rate calculated from eq. 25 assumes that all the soot deposited 
on the surface remains there. In order to determine the actual deposition, the effect of reentrain- 
ment of deposited soot into the gas stream must be considered. The net deposition rate is 
obtained by multiplying the gross deposition rate by a "sticking fraction" that varies between 0 and 
1 . No theoretical models of soot particle adhesion strength were found in the literature. Some of 
the forces responsible for the adhesion of soot to a solid surface include contact potential, electro- 
static, van der Waal's and capillary forces [11]. In the absence of a theoretical treatment of the soot 
adhesion strength, the "sticking fraction" is based on assuming a known statistical distribution for 
the adhesion strength of a soot particle to a wall. A log normal type distribution based on the dust 
particle adhesion work of Berengoltz [12] has been found to also characterize soot particle adhe- 
sion by Makel and Kennedy [30]. In addition, the primary force causing reentrainment and 
reducing deposition has been determined to be the fluid shear stress at the wall [30,31]. Therefore, 
the "sticking fraction" of soot particles is the percentage of particles with adhesion strengths greater 
than the mean fluid shear stress (or alternatively the probability that a given soot particle has an 
adhesion strength greater than the mean shear stress at the wall). 

The adhesion force, F ad , for particles is assumed to be proportional to particle 
diameter, d p , as in eq. 26 

F^kd 2 (26) 

and distributed according to a log-normal distribution in k. The characteristic force is F$o% which 
corresponds to the force required to remove 50% of the deposited particles. The "sticking fraction" 
is defined as one minus the probability that the surface shear stress, x, exceeds the adhesion force. 
Fad. as given in eq. 27. 
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4.2, Model Formulation (cont.) 


R = P(r > FJ) = 



where: k* = xld ^ 

m = ln(K S0%) 
o = standard deviation 


( 27 ) 


The net particle deposition can then be expressed as eq. 28. 

; --<1 -*»./> V (28, 

The Task V test data can not be used to verify the deposition model. The tur- 
bine simulator flow blockage tests yield the percent reduction of C^A (product of the discharge 
coefficient and available flow area.) In the absence of post-test physical measurements of the flow 
area, the only way to relate the pressure drop data to deposit thickness is to assume that the dis- 
charge coefficient is constant. Inspection of post test hardware photographs indicates dramatic 
changes in the shape and roughness of flow passages, indicating that very large variations in Qj 
are to be expected. In addition, the deposition model requires input of the free stream soot volume 
fraction and the adhesion strength distribution. 

Figures 10 and 1 1 show the parametric behavior of the deposition model at 
representative gas generator conditions. The conditions chosen were: 

P c = 1000 psia 

T b = 1520 F 

<j> =0.0001, 0.00005, 0.000025 

v 

The figures indicate two major flow effects on soot deposition. Figure 10 shows the increase in 
deposition mass flux as the flow velocity is increased in various size nozzles. The calculations 
assume a perfectly adherent wall (i.e. R=0). The deposition is proportional to the free stream soot 
volume fraction, <j) e , and approximately linear with velocity (the square of the nozzle diameter). 
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Figure 11. Predicted Effect of Wall Temperature on Carbon Deposition 
for Gas Temperature of 1520 F in a 20:1 Contration Ratio 
Nozzle with L02/RP-1 at MR=0.4 and 2 Ibm/sec. 


4.2, Model Formulation (cont.) 

Figure 1 1 shows the effect of thermophoresis on deposition. As the wall temperature is increased 
towards the free stream temperature the deposition rate can vary by two orders of magnitude. 


4.2.5 1-D Flow Model 


The 1-D flow model is used (1) to convert the time variable in the soot volume 
fraction equation (eq. 14) to an axial position, (2) to calculate the axial presssure profile, and (3) to 
include the effect of various hardware geometries. The model assumes a 1-D isentropic flow, 
with temperature dependent properties and an multiconstituent flow. The equations for this type of 
model are well known [32] and are repeated here for completeness. The principle relations are 
mass continuity, conservation of momentum, and conservation of energy as listed in eqs. 29, 30, 
and 3 1 . 


d , 

— (p, va) = £O i r a 
dx 


(29) 


where: p, = density of ith species 

0} = species production rate 
a = cross sectional area 
r* = throat radius 
x = normalized axial distance 
v = axial velocity 


P 

h 


dv dP 
v — + — 


= 0 


dx dx 

+ -v 2 = H c 
2 


(30) 

(31) 


Using constituent relations for the system enthalpy, pressure, and multi-species gas constant (eqs. 
32, 33, and 34) 


f 


k= I C pi dT + h u 


(32) 
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4.2, Model Formulation (cont.) 

P. = PA t 


( 33 ) 


i=i 


(34) 


the equations can be derived in terms of the gas generator geometry area variation, where the axial 
position and area have been normalized by the throat radius, as listed in eqs 35 through 40. 


dc L = 
dx 

dv _ 
dx 

dr_ _ 
dx 

rfT _ 
dx 
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T 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


Once the bulk flow has been calculated the axial wall shear stress can be 
approximated by existing correlations based on the flow Reynolds number [33] or by 2-D bound- 
ary layer calculation using JANNAF boundary layer codes such as BLIMP-J [34] and Boundary 
Layer Module[32]. 

4.2.6 Combustion Kinetics Model 

In order to calculate the mean axial temperature profile in the gas generator, 
either test data or a chemical kinetics model for the low MR combustion must be used. Simplified 
combustion kinetics models for LO2/RP- 1 such as the Fuel Rich Combustion Model [13] have 
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4.2, Model Formulation (cont.) 

proven adequate for predicting combustion temperature but have not successfully predicted soot 
particle inception and coagulation. For propane, several potential reaction schemes exist which 
should be useful in constructing a kinetics reaction scheme [35]. The completed version of the 
model will use a kinetics model such as the FRCM or a variant. The reaction rates will be 
expressed in terms of modified Amenhius rate expression for j reactions such as eq. 41 . 

k j =ajT n e(~ EalRT ) (41) 

The resulting set of equations are a coupled set of nonlinear ordinary differential equations which 
can be integrated in time with the 1-D flow equations to provide the energy release and temperature 
profile in the gas generator. 

The FRCM starts by assuming that all the oxidizer reacts with a portion of the 
fuel to produce equilibrium combustion products. The JANNAF One-Dimensional Equilibrium 
(ODE) computer program is used to determine the combustion products. The remainder of the fuel 
is presumed to vaporize and react kinetically with the ODE predicted combustion products. The 
JANNAF One-Dimensional Kinetics (ODK) computer program is used for the kinetics calculation. 
The fuel vaporization is modeled using a temperature dependent rate expression to permit its use in 
the the ODK program. The vaporized fuel then undergoes kinetically limited combustion to fomi 
the final products. While this approach is not a true representation of the combustion presumed 
occurring in the gas generator, sufficient adjustment of rate constants and activation energies has 
yielded a reaction set which shows good agreement with RP-1 test data [1]. 

4.3 INJECTOR MIXING EXPERIMENTS 

Cold flow tests were conducted to assess whether typical gas generator combustion is 
diffusion flame-like(mixing limited) or premixed (kinetically limited) relative to typical soot forma- 
tion rates. The tests were conducted using the high injection density injector (Figure 12). Air was 
used to simulate the oxygen and Freon 1 16 was used to simulate the fuel. The tests were 
conducted at the Combustion Research Laboratory at the University of California at Davis. The 
mixing was visualized by imaging the laser light scattered off fuel and oxidizer molecules. The 
Rayleigh scattered light was from a laser sheet passed through the flow and perpendicular to the 
injector face. The propellant simulants were introduced as gases with mixture ratios of 0.25 to 0.5 
and injection momentum ratios comparable to the Task V test conditions. Because the propellants 
were injected as gases, mixing occurs faster than with cold (i.e. below the critical temperature) 
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Figure 12. High Injection Density Injector 
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4.3, Injector Mixing Experiments (cont.) 

supercritical fluid or liquid injection. The flow Reynolds number based on either the injector 
orifice diameter or the chamber diameter was approximately 10 to 100 times lower than that for an 
actual gas generator (which are approximate 10 6 to 10 7 ). but still was fully turbulent. The effect 
of Reynolds number is considered negligible on the tests conducted. 

The measured mixing data (i.e. the relative concentration distributions) indicate large 
regions of unmixed fuel and oxidizer. The observed mixing is less than that required to suppress 
soot particle formation. The conclusion drawn from the test data is that finite rate mixing effects 
on soot particle inception exist with conventional gas generator design. 

4.3.1 Test Apparatus 

The test apparatus is shown in Figure 13 and a schematic diagram is shown in 
Figure 14. The high injection density injector S/N 1202756 was mounted vertically on a test stand 
with the propellant simulants supplied from high pressure supply tanks. The flowrates were set 
using calibrated sonic orifices. The green light output (514 nm) from a 5 watt argon ion laser was 
focused into a sheet 1 .5 in (3.7 cm) wide directed across the flow. A diagnostic chamber (not 
shown in the schematic) was mated to the injector to prevent entrainment of surrounding air. The 
design of the diagnostic chamber is shown in Figure 15. The 12-inch long chamber has two long 
axial slots on opposite sides to permit the sheet of laser light to be introduced across the flow any- 
where along the length of the chamber. At 90 degrees to the sheet are 1.0-inch diameter windows 
located to collect the scattered laser light. The chamber contains four axial passages for heated 
water to maintain a constant chamber temperature (approx. 50 Q during testing to prevent conden- 
sation on the inside windows. 

The laser light scattered off the Freon molecules was imaged onto a gated 
image intensifier with a 2: 1 demagnification ratio. The image intensifier was gated on for 
approximately 10 microseconds using a TTL triggering circuit. The image intensifier provides a 
gain of approximately 1000 to the scattered light intensity. The 10 microsecond “exposure time” is 
sufficiently short to freeze the flow field and provide an instantaneous image of the mixing. The 
output of the image intensifier was imaged onto an RS-170 format video camera which was 
triggered to take a single frame when the intensifier was gated on. The video image was then 
digitized by a VS- 100 frame grabber and then uploaded to an Apollo DN10000 workstation for 
processing and storage. 
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Figure 13. Test Apparatus For Injector Mixing Tests 
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Figure 14. Experimental Apparatus For 2-D Rayleigh Scatter Imaging of Injector Mixing 
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Figure 15. Diagnostic Chamber Assembly (page 1 of 2) 
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4.3, Injector Mixing Experiments (cont.) 


4.3.2 T pst Summary 

A test summary of the injector mixing experiments is shown in Table II. The 
tests simulate mixture ratios from .25 to 0.5. The first five tests used 2-D imaging of the Rayleigh 
scattering from Freon 1 16 to obtain instantaneous maps of the species concentrations. The spatial 
gradients in the concentrations have been used to estimate the mixing. During the second five tests 
the flow was seeded with alumina particles (approximately 1 micron diameter) and laser doppler 
velocimetry measurements of the axial velocity and turbulence intensity at locations of high 
concentration gradient were made. For the LDV measurements the chamber was rotated 90 
degrees and the LDV beams were introduced through the circular window and the doppler signal 
collected in the forward scattering angle by a photomultiplier. The LDV measurements were 
processed using a TSI Burst Doppler Counter and data reduction software. 

Rayleigh scattering of light is a linear process where the scattered light is at the 
same frequency as the incident light. In addition, the intensity of the scattered light is linearly pro- 
portional to the concentration of the scattering molecules in the probe volume (or sheet.) Freon 
1 16 and air were chosen as the simulants because the Rayleigh scattering from the Freon is 
approximately 10 times that for air. The relative concentrations of the fuel simulant were obtained 

using the procedure described below. 

Calibration of Video System 

The variations in pixel sensitivity of the camera and the optical efficiency of the 
lens/intensifier system were calibrated using an image taken with stagnant air. Because the 
concentration is uniform, variations in pixel sensitivity , optics losses, and power variations in the 
laser sheet can be detected on the image. These variations were found to be no greater than 25 %. 
The test images acquired were corrected for this variation by normalizing the data, by dividing the 
digitized intensity reading by the value of the calibration image. 

Determination of Mixture Fraction 

In order to compare the results of the cold flow mixing with previous work on 
soot particle suppression, the relative concentration of the fuel simulant must be expressed in terms 
of a mixture fraction, £ As defined in eq. 42 below, the mixture faction varies between 1 (pure 
fuel) to 0 (pure oxidizer). 
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4.3, Injector Mixing Experiments (cont.) 


{= Z fi IZ f 


z r overall mass fraction of fuel 

= ((1/MR+ DMR)' 1 

Zj.= local mass fraction of fuel 

= % of Full Scale light detected 


( 42 ) 


False color images of the relative fuel concentration at MR=0.3 and MR=0.5 
are shown in Figure 16. The impinging jets of fuel can be seen on each side of the injector. The 
mixing shown is in the axial direction, perpendicular to the injector face. The effect of the injectors 
swirl design can be seen in the upper instantaneous image by the appearance of larger fuel jets on 
the left side of the the images where the fuel is being sprayed into the laser sheet. On the right 
side of the image the fuel jets appear smaller, even though they are the same size as those on the 
left side, because the fuel is being sprayed away from the laser sheet. A large recirculation zone in 
the center of the 'injector is readily apparent. The size of this zone may be exaggerated due to the 
absence of an igniter flow. In general, significant mixing does not appear to occur until at least 20 
injection orifice diameters (approximately 1.5 in) downstream of the injector face. 


Determination of Mixing Rate 


The relevant turbulent mixing parameter for soot formation is the mean scalar 
dissipation rate (the scalar is the mixture fraction £ ) defined in eq. 43, 

x = 2 D(V£) 2 (43) 


where D is the binary diffusion coefficient for air and Freon 116. The scalar 
dissipation rate is a measure of the concentration gradient between the fuel and the oxidizer. A high 
value of the rate implies rapid mixing. The dissipation rate varies throughout the flow field. To 
characterize the entire flow properly by a single average value, a volume weighted average value 
must be used. Such an average can be obtained by analyzing several transverse and axial images 
through the flow. A complete survey of the flow such as this was beyond the scope of this 
investigation. 


An estimate of the average dissipation rate in the primary mixing zone, 
approximately 5 orifice diameters from the injector has been made. The calculated scalar 
dissipation rate at this location is reported in Table II. While no data are available for the 
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Figure 16. 2-D Images of Injector Mixing at MR=0.3 and MR=0.5 


4.3, Injector Mixing Experiments (cont.) 

dissipation rate which causes suppression of soot particles from oxygen/RP-1, Kennedy [24] has 
measured a value of 0.85 sec- 1 for air/propane. The values measured during the cold flow testing 
are approximately 15 to 20% of this value. In addition, LDV measurements of the axial velocity 
were made at the same location as the calculated scalar dissipation rates and they indicate high 
turbulence intensities based on the axial velocity component. The average turbulence intensity, 
which ranged from 35 to 63 % (based on 1200 measurements) are reported in Table II. Typically, 
high turbulence levels correlate with increased mixing. Further investigation of injector designs to 
generate even higher turbulence levels and greater mixing is indicated by the results of these tests. 
Injector designs which substantially promote mixing and short combustion lengths such as a 
HIPERTHIN or swirl vortex designs should provide conditions for less soot formation. 

4.4 MODEL VALIDATION PROGRAM 

A model validation program has been defined (1) to provide data to validate the current 
Low Mixture Ratio Combustion and Carbon Deposition Model, (2) to use the model to design an 
advanced L02/hydrocarbon gas generator, and (3) to use the test data to further the development of 
a CFD based soot formation and deposition model. The time phasing and synergism between 
specific tasks is shown in Figure 17. The empirical variables required to complete and validate the 
current model are shown in Figure 18. The primary data required are (1) axial profiles of the soot 
concentration (i.e. soot volume fraction), (2) soot particle number densities under well-defined 
mixing conditions and (3) soot deposition and soot layer adhesion strength measurements. The 
recommended test program consists of three tiers of testing starting with bench scale testing, then 
high pressure single element testing, and then subscale (1/2 to 1/10 scale) hot fire testing. 

4.4.1 Validation Program Tasks 

Controlled Low Mixture Ratio Combustion Tests 

Currently, little fundamental combustion data exists for very fuel rich oxy- 
gen/RP- 1 and oxygen/propane flames. In particular, the effects of fuel composition (purity and 
additives), fuel pyrolysis reaction paths, soot surface growth rates, and soot deposit characteristic 
are unknown. The initial task of the validation program should be a series of controlled low MR 
combustion experiments to address these issues. The objective of these tests is to obtain sufficient 
data with RP-1 and propane at low mixture ratio conditions so the basic hypotheses and theories of 
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Figure 17. Model Validation Program 
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Figure 18. Data Requirements tor Model Validation 



4.4, Model Validation Program (cont.) 

soot formation and deposition used to develop the model can be verified. The fundamental ques- 
tions to be addressed in the tests are: 

1 . Is the amount of solid carbon formed by liquid fuel coking negligible 
compared to the soot produced by gas phase reactions? 

2. Do the soot particles stop growing at a size of 200 nm, so that ther- 
mophoresis governs the deposition? Do the particles undergo additional 
allogomeration and form particles larger than 1 micron in diameter? 

3. What is the soot surface growth rate for RP-1 and are values which have 
been measured for air/propane at equivalence ratios from 0=1 to 2, valid 
for oxygen/propane at 0= 9 to 12? 

4 . What is the effect of fuel purity on the soot surface growth rate? 

5. What are the major intermediate species and final combustion products 
resulting from low MR RP-1 and propane combustion? 

6. What is the soot layer adhesion strength for RP-1 and propane? 

Insight into these fundamental questions can be achieved through three types of controlled bench 
scale tests. It is recommended that the tests be conducted with RP-1, propellant grade propane, 
research grade dodecane, and high purity grade propane. The latter two fuels should be included 
in the testing to assess fuel blend effects. 

The model uncertainties regarding coking, particle size, and combustion 
species can be reduced by making measurements in the first type of test setup illustrated in Figure 
19, an under ventilated diffusion flame. The required test apparatus consists of an enclosed vessel 
with optical access such as the diagnostic chamber used in the cold flow experiments. A simple 
spray atomizer or coax injection port is used to form a turbulent diffusion flame within the vessel at 
the desired mixture ratio. Because the flame is fuel rich, the flame height is relatively short com- 
pared to a normal diffusion flame and the flame will spread out to the side walls of the chamber as 
illustrated in Figure 19. Optical measurements of soot particle size and gas sampling mea- 
surements can be made throughout the flow. The particle size measurements (Rayleigh and Mie 
scattering) indicate whether the solid carbon is due to soot formation or liquid fuel coking. 
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Figure 19. Controlled Low Mixture Combustion Tests 



4.4, Model Validation Program (cont.) 

The gas sampling measurements using a gas chromatography analysis can identify the species to be 
included in the combustion kinetics model. 

The soot surface growth rate, soot particle inception rate, and fuel composition 
effects can be determined using a counter flow diffusion flame [24,27] as shown in Figure 19. In 
this configuration, fuel and oxidizer flows are directed at each other forming a stationary flame 
front at the stagnation point. The rate of mixing and soot particle residence time can be controlled 
by adjusting the propellant flow rates. The soot surface growth rate can be found by measuring the 
soot volume fraction, particle residence time, and flame temperature. 

The soot deposit adhesion strength can be determined by measuring the soot 
deposition from a hot sooty flow across a cooled cylinder, as shown in Figure 19. The wall shear 
stress can be controlled by adjustment of the flow rate of the hot gas. The cylinder wall temperature 
can be varied by adjusting the cooling on the inside of the cylinder. Measurement of the airborne 
soot volume fraction and the deposition rate can be made using the two color laser measurement 
techniques described in section 4.4.2. From this type of test, both the thermophoretic correlation 
for soot deposition and the soot deposition adhesion strength model can be validated. 

The results from the controlled low mixture ratio combustion tests should pro- 
vide sufficient data to insure that all important processes have been included in the model. In addi- 
tion, some empirical parameters in the model can be quantified. The results should provide suffi- 
cient insight into the remaining uncertainties in the model, such that an efficient hot fire test pro- 
gram can be designed to "fill in" the voids and validate the model. In addition, the data obtained 
should be very helpful in developing a CFD base soot formation and deposition model as described 
below. 


Hot Fire Tests For Model Validation 

While the controlled low MR tests described above are adequate to verify the 
physics of the model and quantify some variables, validation requires direct comparison of the 
model prediction to subscale test data. Two types of tests are recommended for model verification. 
The first type of tests are single element injector tests to screen candidate injector designs and to 
determine the correlation between the mixing, soot number density, and local temperature varia- 
tions. These tests would use a small high pressure single element injector and chamber designed to 
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4.4, Model Validation Program (cont.) 

accommodate optical access such as that used on the Photographic Combustion Characterization of 
LOX/Hydrocarbon Type Propellants, NAS9- 15724 [22] and shown in Figure 20. Single element 
injector testing at realistic injection densities and pressures (1000 to 2000 psia) offers a good com- 
promise between realistic combustion conditions, diagnostic accessibility, and cost. The primary 
objective of these tests should be to develop a database on injector design effects. 

The second type of testing is subscale (1/10 to 1/2 scale) gas generator testing 
such as that which was conducted during Task V of the program. The modular hardware designs 
developed on this program are ideal for the subscale testing. Additional components for making 
the required diagnostic measurement could be incorporated in the test setup. A test matrix would 
be run where the mixture ratio, pressure, and total flow rate would be varied for at least two injec- 
tor designs selected from the single element injector testing. Measurements of the soot volume 
fraction, soot particle size, gas temperature, and deposition at different axial locations would be 
made. The data would provide the basis for anchoring the model. 

Anchor Model and Develop Computer Code 

The results of tests described above would be used to produce a final anchored 
version of the model. The equations for the soot formation, deposition, flow, and combustion 
kinetics would be programmed into a FORTRAN computer code and a detailed user's manual 
would be prepared. Parametric analyses using the model would be conducted to identify gas 
generator designs offering the potential to minimize carbon deposition. 

Design. Fabricate and Test an Advanced LOVHvdrocarbon Gas Generator 

Based on the results of the parametric analysis, the detailed design of an 
advance gas generator should be initiated. The objective of this effort would be to demonstrate the 
use of the model for developing new gas generator designs. Currently, the design of a new gas 
generator is a significant risk because of the uncertainty in extrapolating from the limited number of 
existing designs. To overcome carbon deposition effects, the design of an advanced gas generator 
may be a significant departure from current designs. Testing of the prototype advanced gas 
generator would be conducted and the results compared to the predictions of the model. At the 
conclusion of this effort the anchored code and prototype hardware would be delivered to NASA 
for use as the baseline design in a future L02/hydrocarbon engine development program. 
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Figure 20. Single Element Injector Testing 



4.4, Model Validation Program (cont.) 

CFD Soot Model Development and Supporting Experiments 

Development of a second generation soot formation and deposition model 
using CFD techniques should be initiated. CFD is emerging as a powerful tool for analysis of tur- 
bulent reacting flows in complex geometries. The extension of Level III type models, reviewed in 
Section 4.1.2, from “classical" flame configurations to rocket engine combustion CFD codes 
appears feasible. The objective of this task would be to modify an existing rocket engine CFD 
code (e.g. Aerojet’s BICOMB, FLUENT, etc.) with a submodel to account for soot formation, 
transport, and deposition. The principle task will be to develop a workable turbulent closure 
scheme for the soot formation source terms in the soot particle conservation equation (eq. 1). 
Existing turbulence modeling schemes such as k-e and algebraic stress models are suitable for 
modeling the remaining terms in the equation. A robust CFD based model would reduce the major 
uncertainties introduced by the empirical correlations used in the 1-D model. 

Data to validate the CFD based model can be obtained from the controlled low 
mixture ratio combustion tests and additional support experiments specifically designed support the 
CFD model. The additional support experiments would consist primarily of mapping the 
combustion flow to obtain both time averaged and instantaneous measurements of the velocity, 
temperature, and the injector spray. Such detailed measurements are not required to validate the 1 - 
D empirical model, but are needed to validate the flow field predictions produced by a CFD based 
model. The CFD model offers greater potential applicability and accuracy than the empirical 1-D 
model. The disadvantage of the CFD approach is that it may not yield a reliable gas generator 
design tool in the near future. Alternatively, the empirical 1-D model developed in this study 
should be readily achieved through the test program described above. 

4.4.2 Diagnostics and Hard ware Requirements 

Validation of model will require measurement of the soot parameters, dis- 
cussed in Section 4.3.2, in the controlled low MR combustion tests and the hot fire validation 
tests. In order to measure the parameters indicated in Figure 18, nonintrusive diagnostics will be 
required. A summary of potential diagnostics is given in Table III and described in detail below. 
The diagnostics have been evaluated based on capability, complexity, and use in a hot fire test pro- 
gram. In addition to the nonintrusive diagnostics listed, secondary measurements of pressure, 
temperature, and species concentrations using conventional instrumentation and gas sampling are 
assumed to be available during testing and are not discussed in detail. 
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4.4, Model Validation Program (cont.) 

T .aser Extinction 


The average soot volume fraction along a line-of-sight can be measured from 
the attenuation of laser power of a beam which passes through a sooty gas flow. The extinction 
(attenuation) of the laser beam power is primarily due to the absorption of light by soot particles. 
Rayleigh scattering of laser light from soot particles typically accounts for less than 1 percent of the 
total power loss and its contribution is neglected in determining the soot volume fraction [371. In 
the dense spray region near the face of the injector, attenuation of laser light due to scattering and 
absorption by untamed fuel droplets can occur. The effect of droplets along the beam path lengt 
can be accounted for by measuring the extinction at two separate frequencies. The ratio of the 
attenuated laser beam power, I, to the incident power, Io, is given by eq. 44. 

I/Io=exp(-k„,L) <44) 

where: k„ t = Extinction Coefficient 


L= Path Length 


The procedures for determining the average soot volume fraction along the line-of-sight 
of the beam is to first determine the extinction coefficient, lw The incident laser power, 1 0 , and 
the path length , L, are assumed known as indicated on Figure 21a. The attenuated power, I, is the 
measured with a photodiode and the extinction coefficient is calculated directly from eq. 44. The 
soot volume fraction is related to the extinction coefficient by eq. 45. 


#- 


6/rIm 


(ft ~ 1) 

|V + 2>J 


(45) 


where: X 0 = wavelength of laser 

n= refractive index of soot 


= 1.54-0.57i 

A large attenuation corresponds to a large extinction coefficient and a large soot volume fraction as 
illustrated in Figure 21a. 
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Figure 21. Fundamentals of Laser Diagnostics for Model Validation 
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Figure 21. Con’t 
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4.4, Model Validation Program (cont.) 

Laser extinction has been used extensively for the past 15 years to study soot 
formation. The technique has also been used on large scale combustion devices such as diesel 
engines |38J and gas turbine engines [26J. The technique requires only limited optical access (5 to 
10 mm diameter windows). The laser source for the technique is a relatively inexpensive, 
continuous-wave (cw) argon ion laser with from 3 to 5 watts of power. 

Ravleieh Scattering 

Rayleigh scattering of laser light is an elastic scattering process, the scattered 
light is at the same frequency as the incident light. Rayleigh scattering occurs from both molecules 
and small particles where the wavelength of light, X, and the particle diameter, d p , are Jtdp/X<l . 
For particles where 50>7tdp/X>l, the elastic scattering is governed by Lorentz-Mie theory. For 
soot particles the elastic scattering is typically assumed to be Rayleigh scattering [33, 39], 
Measurement of the Rayleigh scattered light in combination with a laser extinction measurement 
can be used to determine the soot particle number density and the average soot particle size. 

A typical setup for a Rayleigh scattering measurement is shown in Figure 21b. 
Scattered light from a small segment of length, /, of the laser beam is collected by a lens and 
focused into a photomultiplier detector. Light scattered off unbumed fuel and oxidizer droplets 
(typically larger that 10 microns) can be distinguished by monitoring the scattered light at a second 
angle (e.g. 45 degrees). Since Rayleigh scattered light (vertically polarized) scatters equally in all 
directions and Mie scattered light (from large particles) does not scatter uniformly, the Rayleigh 
signal can be extracted by monitoring both detectors and rejecting data when the signal differ. The 
intensity of the Rayleigh scattered light, I, is proportional to the number density of soot particles in 
the probed volume (typically 1 mm3) a s given by eq. 46. 

I = C 1 0 NQ / (ds/dQ) (46) 

where: 

C= Collection Efficiency 
Io=Incident Laser Intensity 
N= Soot Number Density 
Q=Solid Collection Angle 
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4.4, Model Validation Program (cont.) 

1= length of probe volume 
(ds/dQ)= Scattering Cross Section of Soot 
The scattering cross section for soot can be derived as eq.47 [37]. 


4 2 


&jt ( 2nX \ . i 
(ds/dQ) = — y— | Irl 


where: 

/b= laser wavelength 
n= refractive index for soot 
djn 2 -l\ 


( 47 ) 


y = - L 


8 In +2. 


The soot extinction coefficient, k^,, determined from a simultaneous laser extinction 
measurement, is related to the average particle diameter through the extinction efficiency, Q, u , as 
given in eq. 48. 


O = ■ 


K 

N 


v<; 


(48) 


The extinction efficiency is the ratio of the extinction cross section to the geometric cross sectional 
of the soot particle. The extinction efficiency is equal to the sum of the absorption efficiency 


area 


and the scattering efficiency as given in eq. 49. 

Qtxx- Qabs + Qsca 


(49) 


For Rayleigh scattering particles with a complex refractive index the absorption and scattering effi- 
ciencies can be shown as given in eqs. 50 and 51 [37 ]. 


Ryl 

<2^ = ji R e0r ) 


Qsca = 4 (ds / dfi) / ( nd p 2 ) 


(50) 

(51) 
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4.4, Model Validation Program (cont.) 

Equations 47, 48, 49, 50, and 51 can be be solved to determine the average soot particle diameter, 
dp, and the soot particle number density, N. A map of the soot number density through the flow 
can be determined by scanning the beam through the flow and adjusting the focus of the detection 
optics. 

Rayleigh scattering has been used extensively in fundamental soot formation 
studies[21,26] and in combustors. This will be the primary diagnostic technique used to determine 
the soot particle number density at inception and to determine the soot surface growth reaction rate. 

leaser Doppler Velocimetrv (LDV) 

Laser Doppler velocimetry (LDV) techniques have been developed for 
measuring mean and instantaneous velocity components and turbulent velocity parameters. 
Commercially available "turn-key" systems are available and have been applied to a wide range of 
practical combustion devices [36]. A typical cross beam LDV configuration is shown in Figure 
20c. For a given velocity component two polarized, cw laser beams of equal strength and fre- 
quency are focused so that they cross at a point in the flow. The intersection of the beams forms a 
control volume typically 1 mm 3 . Within the control volume, the beams form a fringe pattern of 
light and dark regions due to constructive and destructive interference of the beams. A burst 
Doppler signal as shown in Figure 21c is detected as a particle passes through the control volume. 
The peaks and valleys correspond to the particle being a light or dark region. The overall shape of 
the signal is due to the Gaussian profile of the laser beams. Signal processing electronics are used 
to eliminate measurements occurring when multiple particles are in the control volume or the seed 
particle size is too large or too small. The flow velocity is proportional to the period of the Doppler 
signal as given by eq. 52. 

U=X/(2Tsin0) ( 52 ) 

where: X= laser wavelength 

T= period of Doppler signal 
0= half angle of intersecting beams 

LDV is a mature laser diagnostic technique. Extensive data-processing soft- 
ware is available enabling a very large number of measurements (on the order of 1000) to be taken 
rapidly (1 to 10 seconds depending on the seeding rate) at a given point in the flow. Two and 
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4.4, Model Validation Program (cont.) 


three dimensional systems configured with fiber optic probes are available, where the signal is 
collected though the same optics which deliver the beams. This "back scatter" configuration is the 
preferred arrangement (compared to the forward scatter arrangement used in the cold flow injector 
mixing experiment) for the hot fire test measurements, because only one optical port is required. 


Laser Induced Fluorescence (LIF) 

Laser induced fluorescence spectroscopy (LIF) is a technique used to measure 
the concentration of important combustion radicals species such as OH and CH. In addition, the 
relative concentration of soot particle precursors (polycyclic aromatic hydrocarbons) can be 
measured. LIF consists of an absorption process followed by an emission process. The laser 
beam frequency (typically a high power, tunable, pulsed laser) is tuned to coincide with an 
absorption wavelength of an electronic transition of a molecule of interest. Energy from the laser is 
absorbed by the molecule and an electron is promoted to an excited state. Spontaneous emission 
and collisions of the excited molecules with other molecules causes the transition of excited elec- 
trons to lower energy states and the emission of photons (fluorescence) at the frequency character- 
istic of the transition. The concentration of the particles is determined by measurement of the 
intensity of the fluorescence at specific wavelengths. 

The typical configuration for a point LIF measurement is shown in Figure 21e. 
The laser is tuned to a particular wavelength, Xq (e.g. 306.2 nm for OH,) and is focused at a point 
in the flow. The fluorescence signal, following the laser pulse (typically 5 to 10 mJ and 10 to 50 
nsec in duration) is collected at a 90 degree angle to the beam. The collected light is then spectrally 
dispersed with a grating spectrometer and focused onto a linear diode array. The output of the 
diode array provides an instantaneous measurement of intensity versus wavelength. The technique 
can be extended to planar (2-D) imaging by using cylindrical lenses to expand the beam into a 
sheet and replacing the spectrometer and linear diode array with an intensified video camera. The 
filter is chosen to isolate the fluorescence at a particular wavelength. The video image then 
provides a map of the spatial distribution of the concentration of the particular species being 
probed. 


LIF has been used for mapping flame location in both sooting flows[39] and 
spray combustion [40]. Fuel rich regions of the flow where thermal decomposition of the fuel 
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forms soot precursors can be identified by fluorescence of the CH radical which is present in these 
regions. Fluorescence of the OH under low mixture ratio combustion conditions indicates the 
location of oxidizer rich portions of the flow (i.e. oxidizer streaks). The use of L1F on the single 
element injector tests can provide a hot fire assessment of injector mixing efficiency. Comparison 
of the results for different injector configurations may lead to the design of an injector element 
which minimizes the soot formation rates. 


Two Color Soot Deposition Measurement 

Simultaneous measurement of soot volume fraction and deposition thickness is 
possible using a two color laser diagnostic technique [30]. The technique has been used to mea- 
sure soot deposition rates from hot gas flows to cooled cylinders in crossflow. The technique uses 
two colinear cw laser beams of different frequency. The technique is illustrated in Figure 21d. 

The beams are superimposed on one another and are aligned so that they pass tangent to a surface 
where the deposition will be measured. As soot builds up on the surface, the beams are partially 
blocked. The beams are separated (either by color or polarity) after they pass through the flow and 
the attenuated power is measured. The soot volume fraction can be determined as is described 
above for laser extinction. The additional attenuation of the signal due to physical blockage of the 
beam, which is proportional to the deposition thickness, is determined by comparing the attenuated 
signal of each laser frequency. This technique provides an on-line measurement of the deposit 
thickness during a test. 

Coherent anti-Stokes Raman Spectros copy (CARS) 

Coherent anti-Stokes Raman spectroscopy (CARS) is a very powerful laser 
diagnostic technique capable of measuring the major species and temperature at a point in sooty 
flows. CARS is based on a non-linear, light wave mixing process capable of both high spatial and 
temporal resolution [41]. The technique is illustrated in Figure 21f. Laser beams at frequencies 
coi and 0 ) 2 , called the pump and Stokes, respectively, with a frequency difference selected for a 
particular molecule to be measured, are "mixed" by focusing and crossing the beams in a specific 
manner to generate a laser-like CARS signal beam. Temperature measurements are derived from 
the spectral distribution of the CARS signal beam. Concentration measurements of the probed 
species are determined from the intensity of the CARS signal. 


RPT/E0068.66 


67 


4.4, Model Validation Program (cont.) 

CARS is not considered a feasible diagnostic for subscale testing, but could be 
used for probing the near injector region in the single element testing because of its ability to work 
in the presence of droplets. CARS is by far the most costly and complex (both in terms of 
implementation and data reduction) of the potential diagnostics to be used. 

Hardware Requirements: 

The major hardware requirement for the validation program is the modification 
of existing program test hardware to accommodate optical diagnostics on subscale tests. The 
modular design and known combustion stability characteristics of this hardware provides an ideal 
apparatus for the subscale experiments. A series of “diagnostic flanges” is recommended to pro- 
vide the optical access required by the diagnostics described above. These flanges would be 
designed to be inserted between the existing hardware components. The flanges would be water 
cooled and would provide a nitrogen purge flow over the windows. Three to five flanges could be 
used during the testing to determine the soot volume fraction, soot particle size, velocity, and soot 
deposition at specific locations in the gas generator. 

Figure 22 illustrates a test hardware configuration using the diagnostic flange 
approach. For the case of multiple soot volume fraction measurements shown, a single laser 
source would be used to create multiple beams to be passed through the flow. The attenuated beam 
strength would be measured with a photodiode detector. Laser beam choppers and lock-in ampli- 
fiers would be used to maximize the signal to noise ratio. The signal from the amplifiers would 
then be recorded by the facility data acquisition system. The recorded voltage is proportional to the 
attenuated laser power. 

The conceptual design of a diagnostic flange for simultaneous measurement of 
the soot volume fraction and velocity is shown in Figure 23. The laser beam enters on the left side 
of the flange and is terminated on a photodiode detector. The measured laser power is used to 
determine the soot volume fraction. A commercially available fiber optic, back scatter LDV probe 
is mounted at the bottom of the flange. The fiber optic probe is used to both deliver the beams and 
to collect the Doppler signal. The LDV measurement will require that the propellants be seeded 
with alumina particles. 

The conceptual design of a diagnostic flange for direct soot deposition measurements is 
shown in Figure 24. The beams for the two color deposition measurement technique described 
above can be combined into one beam (with each color light having a different polarization). The 
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Figure 22. Subscale Hardware Configuration for Validation Testing 



Figure 23. Conceptual Design of Diagnostic Flange for Simultaneous Extinction 
and LDV Measurents 


FROM LASER SOURCES 
488 nm & 632 nm LIGHT 



Figure 24. Diagnostic Flange for Soot Deposition Measurement 




4.4, Model Validation Program (cont.) 

single beam is then passed through small (0.125 inch dia.) windows, embedded in a cooled flange, 
on opposite sides of the gas generator. For deposition tests, the beam would pass tangent to the 
leading edge of a cooled surface such as the turbine simulator used in the Task V testing or a 
cooled cylinder in crossflow. The beams would be separated after leaving the gas generator and 
the attenuation for each color beam would be measured. The ratio of the attenuated signals is 
proportional to the soot deposit thickness. The maximum thickness which can be measured can be 
set by using cylindrical lenses to expand the beams into sheets from 5 to 15 mm thick or larger. 
The upper limit on the thickness is limited only by the size of the window access. 
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